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ABSTRACT 

We present MIPS 24/im observations of the Hubble Deep Field South taken with the Spitzer Space 
Telescope. The resulting image is 254 arcmin 2 in size and has a sensitivity ranging between ^12 to 
^30 /xJy rms, with a median sensitivity of ^20 /xJy rms. A total of 495 sources have been cataloged 
with a signal-to-noise ratio greater than 5a. The source catalog is presented as well as source counts 
which have been corrected for completeness and flux boosting. The IR sources are then combined with 
MUSYC optical/NIR and ATHDFS radio observations to obtain redshifts and radio flux densities of 
the sample. We use the IR/radio flux density ratio {qu) to explore the IR-radio correlation for this 
IR sample and find (J24 = 0.71 ± 0.31 for sources detected in both IR and radio. The results are 
extended by stacking IR sources not detected in the radio observations and we derive an average 524 
for redshift bins between < z < 2.5. We find the high redshift (z > 1) sources have an average 524 
ratio which is better fit by local LIRG SEDs rather than local ULIRG SEDs, indicating that high 
redshift ULIRGs differ in their IR/radio properties. So ULIRGs at high redshift have SEDs different 
from those found locally. Infrared faint radio sources are examined, and while nine radio sources do 
not have a MIPS detection and are therefore radio-loud AGN, only one radio source has an extreme 
IRAC 3.6/xm to radio flux density ratio indicating it is a radio- loud AGN at z > 1. 
Subject headings: galaxies: formation — galaxies: evolution — galaxies: starburst 



1. INTRODUCTION 

The tight correlation between the far-IR and radio 
emission of starforming galaxies has been well stud- 
ied in the local universe and establis hed over many or- 
ders of magnitu de in luminosity (e.g. iHelou et al.lll985l : 
lYun et al.lf200l . The radio emission of normal galax- 
ies is dominated by synchrotron radiation from rela- 
tivistic electrons accelerated by supernovae and free-free 
emission from HII regions (ICondod [19921) . Both su- 
pernovae and HII regions are related to young, mas- 
sive stars (M > 10M Q ), and hence radio emission in 
normal galaxies is correlated with star formation. The 
UV emission from young stars is re-radiated by dust 
in the mid-to-far-infrared so IR emission is also related 
to star forming processes. The far-IR/radio correla- 
tion is important because it is a probe of the starform- 
ing processes taking place. Deviations from the cor- 
relation may occur due to, for example, the varying 
strength of galactic magnetic fields, evolution of dust and 



metallicity properties, and increasing Inverse Compton 
losses off the cosmic microwave background fe.g lMurphyl 
120091 lLacki fc Thompson! [20101 ) . The far-IR/radio cor- 
relation can also b e used to select radio loud AGN (e.g. 
iDonlev et al.l [2005 ^) . Thus, it is important to determine 
whether the correlation holds at high redshift. 

A sensitive infrared probe of high-redshift star 
formation is the 24 /im band of the Multiband 
Imaging Photome ter for Spitzer (MIPS) instrument 
dRieke et al.l l2004f) on bo ard the Spitzer Space Tele- 
scope ( Werner et al J [20041 ). Several deep surveys have 
correlated 24 /xm with radio data to study the mid- 
IR/radio relation dBovle et al.ll2007t iBeswick et~aT1 l2008t 
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Ibar et a 1120081: lG arn fc Alexande r 2009). Some of these 
authors (jBovle et al.l [2007: Bcsw ick et al.1 [20081 ) observe 
a change in the mid- IR/radio ratio at faint radio and/or 
24 /im flux densities, but the differences in the results 
are likely to be because of selection biases which lead to 
different populations being studied. 
More recent work in the COSMOS field has found lit- 
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tie variation in the IR-radio ratio of starforming galaxies 
for z < 1.4, and for the limited sample of sources con- 
firmed at z > 2.5 the averag e IR-radio ratio is s till the 
same as that found locally ([Sargent et al.l[20To| ). Dis- 
crepant measurements of the average IR-radio ratio can 
be reconciled by the fact that previous studies selected in 
either IR or radio, o r required a detection in both bands 
(jSargent et al.ll2010j ). 

In 1998 the Hubble Space Telescope (HST) observed 
a region in the southern continuous viewing zone as a 
complementary observation to the Hubble Deep Field 
North (HDF-N. IWilliams et alJll996T >. T his field, known 
as th e Hubble Deep Field South (HDF-S; IWilliams et~aT1 
2000), reaches a limit of ~30 magnitude in the four HST 
WFPC2 broadband filters. Simultaneous observations 
with the other HST instruments, NICMOS and STIS, 
reach similarly deep levels in the near-infrared and ultra- 
violet. In addition to the main HST deep fields, a mo- 
saic of nine flanking fields were imag ed to shallower lev - 
cls of ~25 magnitude in I (F814W) (|Lucas et al.ll2003D . 
Wide-field ground b ased imaging has been performed by 
iPalunas et al.l (|2000H and th e Multiwavelength Su rvey by 
Yale-Chile (MUSYC) team (jGawiser et al.ll200l . reach- 
ing R band magnitude limits of ~25 and 26, respectively. 
This field has been the target of deep radio observa- 
tions reaching ^lOjuJy rm s at 1.4, 2.5, 5.2 and 8.7 GHz 
(|Huvnh et al] 120051 [2007h . 

Confirming the COSMOS results requires a deep and 
reasonably wide MIPS field. In this paper we present the 
MIPS 24/im observations of the Hubble Deep Field South 
and take advantage of the enormous amount of multi- 
wavelength data available to study the IR-radio correla- 
tion of galaxies detected in this field. We use a stacking 
technique in the radio to determine the average IR-radio 
ratio of a flux-limited IR selected galaxy sample. 

In Section 1 we describe the MIPS observations and 
data reduction. We summarize the multiwavelength data 
used for this work, in particular the MUSYC data, in 
Section 2. The MIPS source detection and extraction is 
discussed in Section 4 and here we present the 24 /jm cat- 
alog, source counts and redshift distributions. Section 5 
reports our work on the IR-radio correlation in the HDFS 
and our identification of infrared-faint radio sources in 
the field. We conclude with a summary in Section 6. 
A Hubble constant of 71 kms _1 Mpc _1 , and a standard 
A-CDM cosmology with ft M = 0.27 and Q A = 0.73 is 
assumed throughout this paper. 

2. MIPS OBSERVATIONS AND DATA REDUCTION 

The MIPS 24 /im observations of the Hubble Deep 
Field South were carried out over 2 campaigns. The first 
was in 2004 as part of the Guaranteed Time Observations 
(PI Fazio) and the second in 2007 Cycle 3 observations 
(PLD30873, PI Labbe ). The GTO observations were 
taken using the MIPS Photometry mode and covered six 
pointings in the HDF-S region. The Cycle 3 observations 
consisted of four pointings which covered the area west 
of the WFPC2 field (Figure 1). The coverage (Figure 
1) is very uneven, with effective on sky integration time 
per pixel ranging from about 4500 seconds in the deepest 
area to 900 seconds in the shallow regions. 

The Basic Calibrated Data (BCDs) were downloaded 
from the Spitzer Science Center (SSC) archive and were 
reduced using version 18.1.5 of the MOsaicker and Point 



source Extractor (MOPEX) software from the SSC. As 
recommended, the first frame from each commanded se- 
quence of exposures (i.e. DCENUM=0) was discarded. 
Mosaicing the BCDs at this point resulted in latent ar- 
tifacts. Flatfields were created using the MOPEX rou- 
tine flatfield.pl, which makes a median flatfield from a 
sequence of BCDs. For the GTO data we found that 
bright sources in the field affected the flatfield gener- 
ation because of insufficient dithering within the indi- 
vidual Astronomical Observing Requests (AORs) in the 
GTO data. Therefore we combined all AORs in the GTO 
data to reduce the effect of bright sources being counted 
in the flatfield statistics. All four GTO AORs in 2004 
were observed within four hours of each other, so combin- 
ing all them to make a master fiat is reasonable because 
the instrument flatfield characteristics do not change over 
such a short timescale. The Cycle 3 observations do not 
have the dithering issue, and flatfields were generated for 
each AOR from simply feeding the BCDs of the individ- 
ual AOR into flatfield.pl. 

After flatfielding the BCDs the MOPEX task over- 
lap. pi was used to level the backgrounds. The BCDs 
were then mosaicked in the usual way with MOPEX us- 
ing the default pixel scale of 2.45" per pixel. The outlier 
rejection, which has a default of 5a, was set to the more 
aggressive value of 3a. The resulting image, with units 
of MJy/sr, is shown in Figure 1. This image is about 254 
arcmin 2 in size. 

3. MULTIWAVELENGTH ANCILLARY DATA 

The field was observed by the main three instruments 
on the Hubble Space Telescope in 1998 to depths ri- 
valing that of origina l Hubble Deep Field in the north 
(|Williams et al.l I2000T ). The HST imaging reaches 10ct 
AB magnitudes of appr oximately 27.7 in I band (F814W) 
(jCasertano et al.l 120001 ) in the deepest region. Including 
the flanking fields, Hubble imaging covers only a small 
~50 arcmin 2 area, however. 

More recently, the Multiwavelength Survey by Yale- 
Chile (MUSYC) observed an approximate 37 x 37 ar- 
cmin region surrounding t he Hubble Deep Fiel d South 
in UBVRIz' optical bands (jGawiser et all 120061 ). These 
observations reach UBVRIz' magnitude limits (5a, AB) 
of 26.0, 26.1, 26.0, 25.8, 24.7 and 23.6, respectively. The 
MUSYC team also obtained NIR imaging in J, H and 
K' filters reaching 5a depths of 22.9, 21.8 and 21.1 Vega 
magnitudes over two 10 x 10 ar cmin fields in the ex- 
tended Hubble Deep Field South (|Quadri et al.ll2007t ). 

Infrared Array Camera (IRAC) data from Spitzer was 
also obtained as part of the GTO (PI Fazio) and 2007 
Cycle 3 observations (PID30873, PI Labbe ). The BCDs 
from the SSC archive had the usual procedures applied 
to remove cosmics rays and artifacts, such as column 
pulldown, muxbleed and maxstripe (see lMarchesini et al.1 
2009 for details). Individual background-subtracted im- 
ages were combined into mosaics covering two 10 x 
10 arcmin regions to limiting AB magnitudes (3cr) of 
24.5, 24.0, 22.2, and 22.1 at 3.6, 4.5, 5.8 and 8.0 (j,m, 
respectively (|Marchesini et aLll2009f l . 

Accurate photometric redshifts were determined with 
EAZY from the MUSY C 14 band optical-IRAC data 
(|Marchesini et alll2009f). EAZY, a template-fi tting code 
sinrUartc L jTYPERZ (|Bolzonella et"aiT l2000l) and BPZ 
(|Bemtezl [2000) , has been found, has been found to 
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provide high quality redshifts over < z < 4 
(|Brammer et al.l r2008t ). despite the fact that it is not 
trained on a spectroscopic sample. EAZY has the option 
of linearly combining templates when fitting sources and 
provides multiple estimates of the photometric redshift, 
including z pea k, which is determined from the full prob- 
ability distribution rather than a straight \ 2 minimiza- 
tion. In the case where there are two widely-separated 
peaks in the redshift probability function, z pea k selects 
the peak with the largest integrated probability. In the 
eCDFS the MUSYC team achieved photometric redshift 
accuracies of Sz/(1 + z) ~ 0.008 with EAZY, but this 
result came from supplementing UBVRIz'JHK optical 
imaging, with Spitzer IRAC and 18 medium band optical 
images from Subaru (Cardamone et al. submitted). The 
photometric redshifts from the 14 band optical-IRAC 
data have an accuracy of Sz/(1 + z) = 0.000 ± 0.033 
and the fraction of c atastrophic outliers is only 4% 
(jMarchesini et al.ll2009D . 

The radio observations of the HDF-S were made be- 
tween 1998 and 2001 with the Australia Telescope Com- 
pact Array (ATCA) using all four available frequency 
bands. Between 100 and 300 hours of observing at each 
band yielded images at 1.4, 2.5, 5.2, and 8.7 GHz with 
maximum sensitivities of ^10 /zJy rms ()Huvnh et al.l 
[20051 f2007h . 

4. MIPS IMAGE ANALYSIS 
4.1. Source Detection and Extraction 

Sources detection and extraction was performed using 
the Astronomical Point Source extraction (APEX) tools 
in the MOPEX software. First, APEX performs a local 
background subtraction by taking a median within a box 
surrounding each pixel. The box size was set to 34 pixels, 
or ~14 times the Full Width Half Maximum (FWHM) 
of the 24//m PSF, and the median calculated after re- 
jecting the 100 highest pixels. APEX detects sources 
when a pixel cluster has values which are above the noise 
threshold, but this requires a good estimate of the noise. 
The MOPEX and APEX software can produce several 
types of noise images. The "std" noise image produced 
at the mosaicking stage provides the empirical scatter of 
each pixel from the repeated observations divided by the 
square-root of the number of observations. This under- 
estimates the true noise because it doesn't account for 
pixel-to-pixel correlated noise or confusion noise. APEX 
generates an output "noise" map by calculating the noise 
within a local box, and here we use the same settings as 
for the background subtraction - 34 pixel rectangular box 
and rejec tion of 100 outlier pixels. Using an approach 
similar to IFraver et al.l (|2009f) , the "std" and "noise" im- 
ages are then combined in quadrature with equal weights 
to obtain a combined uncertainty image for source detec- 
tion and extraction. 

We then checked whether the combined uncertainty 
image was reasonable by checking this against an inde- 
pendent total noise value. An independent estimate of 
the total noise in the image was made by measuring the 
flux density in apertures randomly placed on the image 
and calculating the standard deviation of the measure- 
ments. Two hundred apertures were placed randomly in 
60 by 60 pixels regions of a residual image with bright 
sources removed. The standard deviation of the flux den- 



sity measurements agreed with the combined uncertainty 
image to within 50%. 

Sources were detected by APEX if they had a peak 
with a signal-to-noise ratio (SNR) greater than 4 and 
were fitted using the Point-source Response Function 
(PRF) image provided by the SSC. The sources were 
then included in the final catalog if they had a fitted 
SNR greater than 5. This relatively high SNR was cho- 
sen to minimize the number of spurious detections. Even 
with this SNR ratio the initial source list (504 total) was 
visually examined for potential spurious sources such as 
those in the Airy rings of bright sources. Only 9 (~2%) 
sources were removed, suggesting that the APEX algo- 
rithm is reliable. The final 24 /im catalog contains 495 
sources. 

The full 5a catalog is presented in Table [T] A descrip- 
tion of the Table is as follows. 

Column (1) — Source ID. 

Column (2) — Right Ascension in J2000. 

Column (3) — Declination in J2000. 

Column (4) — 24 /im flux density, in fijy. 

Column (5) — 24 /xm flux density uncertainty, in fiiy. 
This is the formal error from APEX. 

Column (6) — Signal to noise ratio of the source. 

4.2. Completeness and Flux Boosting 

The completeness of the source catalog was estimated 
by simulations. Monte-Carlo-like simulations were per- 
formed by injecting 10000 sources at random locations 
on to the map and extracting them using the same tech- 
nique as adopted for the production of the catalog. While 
completeness levels vary significantly across the image as 
a function of coverage and flux density, we can recover 
the overall completeness level of the generated catalog by 
injecting sources over the full image. A single source was 
injected per simulation, and the input flux density varied 
from 50 to 3000 /iJy to sample the full range of interest. 
The completeness as a function of flux density is shown 
in Figure 2. The completeness rises steeply from about 
10% at 100 ^Jy to approximately 90% at 250 /LtJy. The 
50% completeness level occurs at approximately 150 /iJy. 

Source detection algorithms that rely on finding a peak 
above a local noise background lead to flux boosting. 
This is because sources that lie on a noise peak have 
their flux density increased and therefore have a higher 
probability of being detected, while sources which lie on 
a noise trough have decreased flux densities and may 
be excluded altogether. This effect is most pronounced 
in the faintest flux density bins, i.e. in sources with a 
SNR close to the limit of the catalog. The degree of flux 
boosting can be estimated by examining the output to 
input flux density of the simulations (Figure 2). We find 
that flux densities are boosted by about 1% to 4% for 
measured flux densities of 100 to 300 ^iJy, on average. 
The flux boosting is negligible for sources with true flux 
densities brighter than 300 /xJy. 

The simulations can also be used to gauge the posi- 
tional accuracy of the catalog by comparing input and 
output positions. The median of the RA and Dec offsets 
are plotted for various input flux density bins in Figure 3. 
The positional accuracy can be estimated from the stan- 
dard deviation in the offsets of each of these bins. We 
find that at the faintest levels (80 ^Jy) the RA and Dec 
uncertainties are approximately 0.8 and 1.3 arcsec, re- 
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spectively. The positional accuracy is much better than 
0.5 arcsec for sources that are 1 mJy and brighter. This 
is the random component only, and doesn't take into ac- 
count any systematic offsets in the pointing. Spitzer has 
an absolute pointing uncertainty of about 1 arcsec. 

5. RESULTS 

5.1. Source Counts 

The source counts were derived from the SNR > 5 
catalog. We correct the source counts for flux boost- 
ing and completeness using the results of Section 4.2. 
Figure 4 shows the Euclidean normalized source counts 
(dN/dS x S 2 - 5 ) and the results are tabulated in Table 
2. The uncertainties listed are derived from Poissonian 
statistics. The counts from this work span the peak 
near 0.2 - 0.3 mJy. We have also plotted the counts 
from five Guaranteed Time Observation (GTO) fields 
dPapovich et all I2004D. th e First Look Survey (FLS) 
field iMarleau et all (2001 and Spi tzer Wide-field In - 
fraRed Extragalactic (SWIRE) fields IShupe etaLI (|2008h 
for comparison. In the fainter bins HDFS counts are 
broadly consistent with that from previous work, show- 
ing the steep decline for S24 < 0.2 mJy. At brighter flux 
densities the counts are consistent with the FLS, within 
the uncertainties, but lower than that from the GTO or 
SWIRE fields. At -0.5 mJy the SWIRE counts are about 
40% higher the HDFS, and about 20% higher than the 
FLS or GTO source counts. This illustrat es that there 
is sign ificant field to field variation, which IShupe et all 
(2008) estimated to be about 10% for sub-mJy flux den- 
sities in regions as large as the SWIRE fields (> 7 de- 
grees). 

Also overlaid in Figu re 4 are the norm alized 
counts from the mo del of lLagache et all (|T004) . Sec 
lLagache et al.l (|2004l ) for details, but briefly, this model 
uses starburst and normal galaxy templates and fits for 
the evolution of the local luminosity function that is best 
constrained by observed source counts and redshift dis- 
tributions. We note that this model used counts at 15, 
60, 170 and 850 urn from literatu re as well as the 24 /mi 
counts from lPapovich et al"1(|2004[ ). The model is broadly 
consistent with the HDFS, FLS and GTO counts and re- 
produces the turnover at —0.1 mJy that is seen in these 
fields. 

5.2. Optical Counterparts 

We used the MUSYC catalogs to identify optical coun- 
terparts to the MIPS sources and obtain redshift in- 
formation from either the spectroscopic or photometric 
catalogs. The MUSYC imaging overlaps almost the en- 
tire MIPS imaging, with only the — 24 arcmin 2 region 
north of Dec = -60.4994 falling outside the MUSYC 
observations. Thus, 25/495 (5%) MIPS sources do not 
have MUSYC coverage. The MUSYC K-selected cat- 
alog (Quadri et al. 2007) and vl.O optical catalog 
(|Gawiser et all I2006T) were searched for counterparts to 
the rest of the MIPS sample. 

With the —1 arcsec seeing of the observations and the 
astrometric uncertainties of the USNO-A2 stars used for 
astrometry, the positional uncertainties of the MUSYC 
sources are about 0.2 - 0.3 arcsec. For MIPS the la 
positional uncertainty can be estimated by 0.5#fwhm/ 
(S/N) assuming Gaussian noise, and hence a beam size 



of 6 arcsec implies a positional uncertainty of 0.6 arcsec 
(in each coordinate) for the faintest sources. Including 
the Spitzer 1 arcsec absolute pointing uncertainty, the 
lcr offset between MIPS and MUSYC sources should be 
about 1.3 arcsec for the faintest MIPS sources. 

Figure 5 shows a plot of distribution of the offsets 
between all MIPS and MUSYC K-selected (Quadri et 
al. 2007) sources within 2 arcsec. As expected, there 
is a peak in the distribution of offsets at <0.5 arcsec. 
The contribution from chance matches due to the opti- 
cal source density is also shown in Figure 5. We find that 
the number of chance matches meets the number of ob- 
served matches at a radius of ~ 1.5 arcsec and therefore 
use this as our matching radius. 369 MIPS sources have 
a counterpart in the K-selected cat a logs a nd a photomet- 
ric redshift from IMarchesini et all ([2009) . A further 37 
have counterparts in the optical MUSYC catalogs, which 
cover a larger area than the K band imaging. For these 37 
sources EAZY photometric redshifts were calculated from 
the MUSYC UBVRIz optical data plus a narrow band 
5000 A image. A comparison of these seven-band pho- 
tometric redshifts with 166 spectroscopic redshifts from 
literature finds only 18% have Sz/(1 + z) > 0.2. Exclud- 
ing these outliers, we find 5z/(l + z) = 0.027 ±0.046, and 
thus conclude the photometric redshifts are statistically 
reliable. In total 406/495 (82%) of the MIPS sources 
have a MUSYC optical or NIR counterpart. Three have 
more than one MUSYC source within 1.5 and in these 
cases the optical candidates were individually examined 
and t he most likely candidate chosen using the P statistic 
(e.g. iDownes et al|[T98llLillv et allll999f ). 

Spectroscopic redshifts from the literature are available 
for 22 of the MIPS sources, and we assume the photomet- 
ric redshift, where available, for the remaining sources. 
The redshift distribution of the MIPS sources peak at 
z — 0.5 with a tail to high redshift (z > 2, Figure 5). The 
total infrared luminosity (8 - 1000 /mi) of the sources 
is estimated by fitting the 24 /mi flux density to the 
lumin osity dependent SED templates of I Chary fc Elbazl 
([20011) . We integrate the best fit template over 8 - 1000 
/mi to derive the total IR luminosity (LIR) (Figure 6). 
We find luminous IR galaxies (LIRGS, 1O U £ < LIR 
< 10 12 L Q ) are detected out to z ~ 1.1, while ultralu- 
minous IR galaxies (ULIRGs, LIR > 10 12 Lq) are de- 
tected out to h i gher r edshifts. We note that while the 
iCharv fc Elbazl ([200 ID templates exhibit observed FIR 
flux density ratios which are more representative of z ~ 1 
galaxies than Dale and Helou (2002) or Lagache et al. 
(2003) templates (Magnelli et al. 2009), the choice of 
templates can lead to ~ 0.3 dex variation in luminosity. 
Fitting to the 24 /mi flux density alone can give a good 
LIR estimate up to about z — 1.4, but for 1.4 < z < 2.5 
the LIR can be overestimated by up to a factor of four 
(Murphy et al. 2009). 

6. DISCUSSION 

6.1. Infrared- Radio Correlation 

Radio counterparts to the MIR sources were ob- 
tained using the 1 .4 GHz radio catalog of the HDFS 
(|Huvnh et al.l 12005). We applied a simple radial match 
of 3 arcsec, and find 84 MIPS sources have a radio coun- 
terpart. The matching distance is large enough for ob- 
taining bona-fide matches since the 24 /mi positional un- 
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certainties are ~ 1.2 arcsec and the radio positions are 
better than about 1 arcsec in most cases (Huynh et al. 
2005). The P-statistic suggests the probability of one or 
more radio sources within 3 arcsec of a IR one by chance 
is ~ 0.3%, but this estimate does not take into account 
the individual positional uncertainties and assumes the 
radio source population is not clustered. 

The IR/radio correlation as a function of redshift is 
presented in Figure 7, shown as the observable parame- 
ter g24 = log(5 , 24 M m)/5'i.4GHz)- For the 84 sources with 
detections in both the IR and radio the median g24 = 
0.71 ± 0.31. This is co nsistent with th e q 2i = 0.84 ± 
0.28 found in the xFLS lAppleton et~aT1 (|2004[ ) and q 2i 
= 0.66 ± 0.39 in the S ubaruXMMNewton Deep Field 
fSXDF. llbareTallllOO^ . Overl aid on Figure 7 are SEP 
tracks from local templates by Ch ary fc Elbazl (|2001f ). 
We find that up to redshift one the IR sources appear 
to follow tracks of local SEDs, suggesting the IR-radio 
correlation holds to this redshift. Above redshift one the 
24 (im band starts to observe the MIR regime (< 12/im 
rcstframe) where PAH emission and silicate absorption 
strengths vary widely in individual galaxies, and hence 
a wide range of MIR/radio ratios is possible. The high 
redshift (z > 1) sources in Figure 7 also fall within the 
ratios allowed by local templates, but most are best de- 
scribed by the SEDs of relatively low luminosity SEDs 
(L < 10 12 Lsun). These sources have luminosities that 
are ULIRG-like (> 10 12 Lsun, Figure 6) but exhibit ra- 
tios similar to local ~ 10 11 Lsun LIRGs, implying that 
the SEDs of high redshift ULIRGs differ to that of lo- 
cal ones. We speculate that these ULIRGs have strong 
PAH emission, as s een in other IR samples of high red- 
shift ULIRGs (e.g. lYan et all 120071) . but IR spectra are 
needed to confirm that hypothesis. 

To extend the results we stack the IR sample in the 
radio data. First, the IR sources were binned into 6 
redshifts bins, as shown in Table 3. For each redshift 
bin, the radio data was stacked at positions of all un- 
detected sources and a noise-weighted mean image was 
produced. Stacks of width 60" wide were produced and 
the photometry performed by fitting a point source to 
the stacked radio image. All stacks, shown as postage 
stamps in Figure 8, have a significant detection. The re- 
sulting ^24 values from the stacked radio flux densities 
are shown as blue squares in Figure 7. They have been 
combined with the actual detections to derive an aver- 
age <;24 for all MIPS sources in the various redshift bins 
(green squares). We find that the average observed q 2 4 
at high redshift (z > 1) is approximately 0.86 (Table 3) 
and this is consistent with local LIR ~ 10 11 LIRG SED 
(Figure 7). At these redshifts our sample have ULIRG 
luminosities or greater, so this is further evidence that 
the average IR luminous galaxy at high redshift is better 
described by a lower luminosity local template. This re- 
sult is consistent with complementary studies which have 
found that some high redshift IR luminous galaxies have 
cooler SEDs than lo c al galaxies of simil ar luminosities 
(jSevmour et al.|[20Tot IMuzzin et al.|[20loh . 

As mentioned before, the high redshift ULIRGs may 
exhibit stronger PAH emission than local galaxies, lead- 
ing to a greater q 2 4 then expected from local SEDs. How- 
ever, the increased q 2 4 at high redshift may also result 
from suppressed radio emission. If a galaxy has just be- 



gun a burst of star formation then the MIR dust emission 
may be enhanced due to star formation but the super- 
nova rate, which powers the radio emission, lags the star 
formation by about ~ 30 Myr (the approximate lifetime 
of 8 M Q stars). It is possible that some of the sources 
at high redshift are nascent starbursts. Increased inverse 
Compton scattering of the cosmic microwave background 
at high redshift can lead to lower radio emission than in 
the local universe, and since the energy density of the 
CMB scales as (1 + z) 4 this could be significant at high 
redshifts. Another potential explanation is that the high 
redshift galaxies have lower magnetic field strengths than 
local galaxies, leading to decreased electron confinement 
and therefore lower radio emission. 

6.2. Infrared- Faint Radio Sources 

Infrared-Faint Radio Sources (IFRSs) are a class of ra- 
dio sources recently discovered in th e Australia T el escop e 
Large Area Survey (ATLAS) by iNorris et all (|2006D . 
These are radio sources brighter than a few hundred /iJy 
at 1.4 GHz which have no observable infrared counterpart 
in the Sp itzer Wide-area Infra red Extragalactic Survey 
(SWIRE: [Lonsdale et al~ll2004[ ). Most have flux densities 
of a few hundred /xJy at 1.4 GHz, but some are as bright 
as 20 mjy, resulting in extreme infrared to radio ra- 
tios. Recent VLBI detections have constrained the radio 
sources sizes to les s than 0.03 arcsec. suggesting IFRSs 
are c ompact AGN (jNorris et al.ll2007t iMiddelberg et al.1 
2008). Deeper Spitzer le gacy survey data in the extended 
Chandra Deep Field South yielded two IRAC detections 
of the four IFRSs in that field, and SED modeling of 
these sources with the new constraints sho wed that they 
are c onsistent with z > 1 radio-loud AGN (jHuvnh et al.l 

l2oTnh . 

The MIPS data of the Hubble Deep Field South al- 
lows us to search for IFRSs in this field. There are 
nine ATHDFS sources with Siaghz > 200 fiJy that 
do not have a MIPS counterpart. These are listed in 
in Table 4 along with their MIPS (5a) limit and IRAC 
properties. All of these sources have low q 2 4 values con- 
sistent with radio-loud AGN, and most have moderate 
photometric redshifts of 0.417 < z < 1.548. Most are 
also relatively bright in the IRAC channels, and would 
have been detected by SWIRE depth IRAC imaging. Of 
the nine sources, one (ATHDFS.J2 23205.9-603857) is 
not in the 5cr IRAC catalog of iMarchesini et al.l (2009) 
and one is not in the IRAC field. Hence only one 
source matches the Norris et al. 2006 criteria for 
IFRS. One IFRS in the 252 arcmin 2 MIPS HDFS re- 
gion is consistent with the IFRS source density of ~ 
16/deg 2 in the eCDFS. ATHDFS_J223205.9-603857has 
a S^.Q^m (5cr) limit of 2 /iJy and hence a radio-to- 
IRAC flux density ratio (Siaghz/ S^.e^m) greater than 
about 130. Similarly high radio-to-lRAG flux density 
ratios are found for both IFRS and high redshift ra- 
dio g alaxies (HzRGs) (|Sevmour et aLll2"o"07t iHuvnh et al.l 
2010), suggesting they are high-redshift AGN. For 
ATHDFS_J223205.9-603857, the low value of g24 (< 
—0.50) and high value of Siaghz/ S^.e^m therefore sug- 
gests it is a high redshift (z > 1) radio-loud AGN. 

The radio spectral index between 1.4 and 2.5 GHz for 
these nine ATHDFS sources is also listed in Table 4. All 
but two have radio SEDs with a spectral index a ~ —0.8 
(S oc v a ), which is consistent with synchrotron emission. 
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The IFRS candidate ATHDFS_J223205.9-603857 has a 
spectral index a = —0.80 ±0.29. This result is consistent 
with Middelberg et al. 2010 (submitted) who find IFRSs 
have a < —0.7. 

7. SUMMARY 

We have presented the MIPS 24 /im observations of 
the Hubble Deep Field South from the Spitzer Space 
Telescope. This new MIR image, 250 arcmin 2 in area, 
reaches a maximum sensitivity of ~12 ^Jy rms, with the 
bulk of the image achieving ~ 20 /iJy rms. A catalog of 
sources brighter than 5cr was compiled, comprising 495 
sources. The completeness corrected and flux-deboosted 
24 /jm source count was presented, and while the counts 
are consistent with results in the literature they show 
that there is significant field to field variation. 

The MIPS sources were matched to MUSYC pho- 
tometric redshift and ATHDFS 1.4 GHz radio cata- 
logs. The IR-radio ratio q 2 4, log(S , 24 A im)/S'i.4GHz), for 
84 sources with both IR and radio detections have a me- 
dian q 2 4 = 0.71 ± 0.31, consistent with previous results. 
We find that while there is diversity in galaxy SEDs, 
the IR-radio correlation appears to hold out to z ~ 1. 
The results were extended by splitting the sample into 
redshift bins and stacking all IR sources not detected 
in the radio to derive a an average radio flux density 
and hence an average q 2 4- The observed average <?24 was 
compared with that expected from local SED templates 
and at z > 1 the observed average 524 is found to be 
better described by LIRG templates, even though these 
sources have ULIRG-likc luminosities. This is further 
evidence that the average IR luminous galaxy at high 
redshift is better described by a lower luminosity local 



template. The change in (724 may be due to enhanced 
MIR emission, possibly from hot dust heated by AGN 
or relatively bright PAH emission, or suppressed radio 
emission due to lower strength magnetic fields, inverse 
Compton scattering off the CMB. Several of the high 
redshift galaxies may be nascent starburst where the ra- 
dio emission from supernovae is delayed compared to the 
IR and hence a high </24 results. More detailed IR ob- 
servations and spectroscopy to determine the AGN/SB 
nature of the IR galaxies is required to differentiate be- 
tween these possible scenarios. 

The Herschel Space Observatory, launched in May 
2009, is currently undertaking deep IR surveys in bands 
ranging from 70 to 450 fim. These bands are close to the 
peak of the IR emission from dust and not significantly 
affected by AGN or PAH emission. With its unprece- 
dented depth and resolution, Herschel will observe large 
complete samples of IR galaxies to redshifts of z = 2 and 
beyond. In combination with upcoming radio facilities 
such as eVLA, ASKAP, and MeerKAT, this will allow a 
major improvement in the study of the IR-radio correla- 
tion in the high redshift universe. 
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JPL/Caltech. This material is also based on work sup- 
ported by the National Science Foundation under grant 
AST-0807570. 



REFERENCES 



Applcton, P. N, et al. 2004, ApJS, 154, 147 
Bem'tez, N. 2000, ApJ, 536, 571 

Beswick, R. J., Muxlow, T. W. B., Thrall, H., Richards, A. M. S., 

& Garrington, S. T. 2008, MNRAS, 385, 1143 
Bolzonella, M., Miralles, J., & Pello, R. 2000, A&A, 363, 476 
Boyle, B. J., Cornwell, T. J., Middelberg, E., Norris, R. P., 

Appleton, P. N., & Smail, I. 2007, MNRAS, 376, 1182 
Brammer, G. B., van Dokkum, P. G., & Coppi, P. 2008, ApJ, 

686, 1503 

Casertano, S., et al. 2000, AJ, 120, 2747 
Chary, R., & Elbaz, D. 2001, ApJ, 556, 562 
Condon, J. J. 1992, ARA&A, 30, 575 

Donley, J. L., Rieke, G. H., Rigby, J. R., & Perez-Gonzalez, P. G. 

2005, ApJ, 634, 169 
Downes, A. J. B., Peacock, J. A., Savage, A., & Carrie, D. R. 

1986, MNRAS, 218, 31 
Frayer, D. T., et al. 2009, AJ, 138, 1261 
Garn, T., & Alexander, P. 2009, MNRAS, 394, 105 
Gawiser, E., et al. 2006, ApJS, 162, 1 

Helou, G., Soifer, B. T., & Rowan-Robinson, M. 1985, ApJ, 298, 
L7 

Huynh, M. T., Jackson, C. A., & Norris, R. P. 2007, AJ, 133, 1331 
Huynh, M. T., Jackson, C. A., Norris, R. P., & Prandoni, I. 2005, 
AJ, 130, 1373 

Huynh, M. T., Norris, R. P., Siana, B., & Middelberg, E. 2010, 

ApJ, 710, 698 
Ibar, E., et al. 2008, MNRAS, 386, 953 
Lacki, B. C, & Thompson, T. A. 2010, ApJ, 717, 196 
Lagache, G., et al. 2004, ApJS, 154, 112 



Lilly, S. J., Ealcs, S. A., Gear, W. K. P., Hammer, F., Le Fevre, 
O., Crampton, D., Bond, J. R., & Dunne, L. 1999, ApJ, 518, 
641 

Lonsdale, C, et al. 2004, ApJS, 154, 54 
Lucas, R. A., et al. 2003, AJ, 125, 398 

Marchesini, D., van Dokkum, P. G., Forstcr Schreiber, N. M., 

Franx, M., Labbe, I., & Wuyts, S. 2009, ApJ, 701, 1765 
Marleau, F. R., et al. 2004, ApJS, 154, 66 

Middelberg, E., Norris, R. P., Tingay, S., Mao, M. Y., Phillips, 

C. J., & Hotan, A. W. 2008, A&A, 491, 435 
Murphy, E. J. 2009, ApJ, 706, 482 

Muzzin, A., van Dokkum, P., Kriek, M., Labbe, I., Cury, I., 

Marchesini, D., & Franx, M. 2010, ArXiv e-prints 
Norris, R. P., et al. 2006, AJ, 132, 2409 

Norris, R. P., Tingay, S., Phillips, C, Middelberg, E., Deller, A., 

& Appleton, P. N. 2007, MNRAS, 378, 1434 
Palunas, P., et al. 2000, ApJ, 541, 61 
Papovich, C, et al. 2004, ApJS, 154, 70 
Quadri, R., et al. 2007, AJ, 134, 1103 
Rieke, G. H., et al. 2004, ApJS, 154, 25 
Sargent, M. T., et al. 2010, ApJS, 186, 341 
Seymour, N., et al. 2007, ApJS, 171, 353 

Seymour, N., Symeonidis, M., Page, M. J., Huynh, M., Dwelly, 

T., McHardy, I. M., & Rieke, G. 2010, MNRAS, 402, 2666 
Shupc, D. L., et al. 2008, AJ, 135, 1050 
Werner, M. W., et al. 2004, ApJS, 154, 1 
Williams, R. E., et al. 2000, AJ, 120, 2735 
Williams, R. E., et al. 1996, AJ, 112, 1335 
Yan, L., et al. 2007, ApJ, 658, 778 

Yun, M. S., Reddy, N. A., & Condon, J. J. 2001, ApJ, 554, 803 



7 



TABLE 1 
The MIPS 24/ai catalog. 



ID 




RA 


Dec 






SNR 


MUSYC.MIPSJ223100.6- 


-603403 


337.7524414 


-60.5674782 


882.0 


46.2 


31.8 


MUSYC.MIPSJ223105.6- 


-603357 


337.7733459 


-60.5657158 


355.8 


15.0 


15.2 


MUSYC.MLPSJ223106.3- 


-603421 


337.7760620 


-60.5725441 


719.7 


15.1 


32.2 


MUSYC.MIPSJ223110.7- 


-603126 


337.7944946 


-60.5239792 


717.2 


55.2 


28.5 


MUSYC.MIPSJ223111.0- 


-603518 


337.7959290 


-60.5884361 


205.6 


11.3 


10.2 



■ Table 1 is available in its entirety via the online version. 



Range in S 2 4 M m (&4^) N corrected N (dN/dS)S 2 - 5 
(/fJy) (AtJy) [mjy 1 - 5 deg~ 2 ] 



110 


- 134 


115 


24 


121.4 


368.8 ± 75.3 


134 


- 163 


145 


65 


149.4 


650.7 ± 80.7 


163 


- 198 


174 


70 


106.8 


585.4 ± 70.0 


198 


- 241 


213 


83 


97.5 


738.0 ± 81.0 


241 


- 293 


254 


70 


76.0 


753.0 ± 90.0 


293 


- 356 


316 


54 


58.0 


788.1 ±107.2 


356 


- 433 


388 


37 


38.4 


706.3 ± 116.1 


433 


- 526 


454 


29 


29.7 


659.0 ± 122.4 


526 


- 640 


590 


13 


13.2 


444.6 ± 123.3 


640 


- 779 


746 


14 


14.3 


712.5 ± 190.4 



TABLE 2 

The MIPS 24/um source count for Hubble Deep Field South. N is the raw number of galaxies in the bin, and 

CORRECTED N IS THE NUMBER CORRECTED FOR INCOMPLETENESS. THE MEDIAN FLUX DENSITY OF THE BIN, (SWpm), HAS 

BEEN CORRECTED FOR FLUX BOOSTING. 



z 


median z 


Ndet 


N s tack 


stacked SiAGHz (juJy) 


stacked 524 


average 524 


< z < 0.25 


0.149 


8 


25 


15.2 ± 4.1 


1 OB+U.1U 


i - uy -0.07 


0.25 < z < 0.50 


0.389 


19 


54 


23.8 ± 3.5 


i 1O+0.06 
U1S-0.07 


n qo+0.05 

u-y^-o.os 


0.50 < 2 < 0.75 


0.580 


14 


65 


22.8 ± 5.2 


1 05+ 009 

i - uo -o.n 


Q1+ 07 


0.75 < z < 1.00 


0.890 


3 


57 


16.4 ± 5.0 


i 14+ - 12 


1 14+O.IO 


1.00 < 2 < 1.50 


1.270 


10 


72 


18.2 ± 3.0 


i -,1+0.07 
1 - 11 -0.08 


n sfi+ ' 06 


1.50 < 2 < 2.50 


1.787 


13 


61 


12.0 ± 2.9 


-1 97 +0.10 
i - z '-0.12 


0-»8±oloe 


no redshift 




16 


73 


18.7 ± 3.6 


, 19 +0.08 
± - lz -0.10 


0.89ir « 



TABLE 3 

Summary of radio stacking results of MIPS sources. N det is the number of MIPS sources in each bin which is 

DETECTED IN THE RADIO IMAGE. N stack IS THE NUMBER OF MIPS SOURCES NOT DETECTED IN THE RADIO IMAGE AND 
THEREFORE THE NUMBER STACKED IN THIS ANALYSIS, RESULTING IN THE STACKED SiAGHz LISTED IN COLUMN 5. STACKED 
§24 IS THE IR-RADIO FLUX DENSITY RATIO FOR THE STACKED SAMPLE. AVERAGE §24 IS THE IR-RADIO FLUX DENSITY RATIO 
FOR ALL MIPS SOURCES IN THE REDSHIFT BIN, OBTAINED BY COMBINING THE DETECTIONS WITH THE STACKED RESULTS. 
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rflydio ID 






0\.&GH z 






g24 




Zphot 


'2.5GHz a 
"l.iGHz — 








(mjy) 


III lv~) 












ATTTDF^ 199^994 ^ 


—604113 


130 


4.108 


112 




— 1.50 


37 


0.746 


— n 48 4- 09 


ATHFJF^ 199^94^ 


— UUOOJ ( 


100 


0.843 


31.2 




—0.93 


27 


1.548 


_n «q -L n Dfi 
— u.uy zn u.uu 


ATTTFlPC; 799^90^ Q 


— ouooo / 


70 


n 9^a 


< z 




— U.OD 


J> lH 


b 


n so -(- n oq 
— u.ou zc u.zy 


ATHDFS.J223350.0 


-603741 


80 


0.670 


156 


< 


-0.93 


4.4 


0.64~ 


-0.36 ± 0.05 


ATHDFS.J223232.8 


-603737 


65 


0.645 


92.6 


< 


-1.00 


7.0 


1.13 


-0.51 ± 0.08 


ATHDFS.J223350.5 


-603503 


110 


1.252 


>4_c 


< 


-1.06 


<310^ 


0.417 


-0.85 ± 0.05 


ATHDFS.J223258.5 


-603346 


95 


1.010 


21.31 


< 


-1.03 


47.4~ 




-0.65 ± 0.05 


ATHDFS.J223308.6 


-603251 


75 


0.821 


101 


< 


-1.04 


8.1 


0.746 


-0.86 ± 0.03 


ATHDFS.J223113.5 


-603147 


115 


0.831 


d 


< 


-0.86 


d 


b 





TABLE 4 

Summary of radio sources that have no MIPS 24 /zm counterpart and hence are candidate infrared-faint 
radio sources. Only ATHDFS_J223205. 9— 603857 is faint enough at 3.6/im to meet the original IFRS definition 
of Norris et al. 2006. All of these sources have 1724 values suggestiv e of radio-loud AGN. The 3.5 ^m data is 

from Marchesi ni et al] (|2009D . 



a Radio spectral indices are from lHuvnh et al.l 120071 ') ■ and S 
b Not detected in MUSYC optical imaging. 

C IRAC image of source has artefacts, so photometry is inaccurate. 
d Out of IRAC coverage. 




338.100 338.000 337.900 337.800 



Fig. 1. — LEFT: The MIPS 24 micron image of the Hubble Deep Field South region. The green chevron figure marks the area covered by 
the original HST WFPC2 deep field, the squares to the east and south are the HST STIS and NICMOS deep fields, respectively. RIGHT: 
The coverage map of the MIPS 24 micron image, white indicates most coverage and black indicates no coverage. 
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Fig. 2. — LEFT: Completeness as a function of input flux density, as derived from the Monte-Carlo simulations. Completeness is the 
number of extracted sources divided by number of input sources. RIGHT: The distribution of (input flux density - output flux density) /input 
flux density as a function of output flux density for the simulated sources. The solid red line is the mean of the simulation and the dashed 
lines mark the 1 sigma upper and lower bounds. 
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Fig. 3. — The offset in RA (left) and Dec (right) between the recovered positions of sources in the simulation and the true input positions, 
as a function of input flux density. The error bars mark the 1 sigma uncertainty in the position as a function of input flux density. 
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Fig. 4. — The Euclidean normalized differential 24 /im source counts in the Hubble Deep Field South (black squares). These counts have 
been corrected for incompleteness and flux boosting as described in Section 4.2. The solid line shows the prediction from the model by 
Lagache et al. 2004. Counts from other deep surveys in the literature are plotted for comparison. 
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Fig. 5. — LEFT: Histogram of the total offsets between the MUSYC optical candidate and the 24 (im source. Most MIR sources have an 
optical counterpart within 1 arcsec. The solid line is the expected number of coincident sources with the same offset, assuming the MUSYC 
source density of 56/arcmin 2 . RIGHT: The redshift distribution of the 406 MIPS sources in HDFS with a MUSYC optical counterpart. 
There are spectroscopic redshifts for 22 sources. 
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Fig. 6.— Total IR luminosities of the 358 MIPS sb&fiesl with rcdshift information. They were derived by fitting to the luminosity 
dependent SEDs of Chary and Elbaz (2001). Inset The corresponding IR luminosity histograms of the sample for z < 0.5 (black), z < 1 
(grey) and z < 2 (white). 
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Fig. 7. — The observed 24 (im versus radio flux density ratio (924) as a function of redshift. HDFS sources detected in MIR and radio are 
marked as black circles. MIR only sources not detected in the radio are plotted as lower limits assuming the conservative radio detection 
limit of 100 fijy. Sources with no redshift information are plotted at z = 3.8. Red circles are results from the xFLS (Frayer et al. 2006). 
The blue squares are average 524 for the MIR sources undetected in the radio, where the average radio flux density is derived from a 
stacking analysis. The green squares are the average 524 for all MIR sources, i.e. combining the blue squares with radio detections. The 
dotted blue lines are SED tracks from Chary and Elbaz (2001) for galaxies with total infrared luminosities of 10 9 , 10 11 , 10 12 and 10 13 Lq, 
going from top to bottom. 




Fig. 8. — Postage stamps of the radio stacks of MIPS sources, for the redshift bins as shown. The stacks are lxl arcmin in size. In all 
cases there is a significant radio flux at the center of the stack, which is reported in Table 3. 



